Thermodynamic geometry of a kagome Ising model in a magnetic field 
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Abstract 

We derived the thermodynamic curvature of the Ising model on a kagome lattice under the presence of an 
external magnetic field. The curvature was found to have a singularity at the critical point. We focused on 
the zero field case to derive thermodynamic curvature and its components near the criticality. According to 
standard scaling, scalar curvature R behaves as |/3 — /3 c | a ~ 2 for a > where /3 is the inverse temperature and 
a is the critical exponent of specific heat. In the model considered here in which a is zero, we found that R 
behaves as |j8 — j8 c | a_1 . 
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1. Introduction 

Statistical thermodynamics provides the basic 
tools needed for understanding macroscopic phe- 
nomena in the physical world. Among the geometric 
representations of statistical mechanics and thermo- 
dynamics of physical systems known are included 
three kinds of Riemannian metric structures to rep- 
resent fluctuation theorem that have been introduced 
by Weinhold and Ruppeiner [QIS]- One of these met- 
rics depends on derivatives of internal energy, while 
the other depends on entropy. Another useful met- 
ric was called Fisher-Rao metric fl. Useful infor- 
mation can be gained about a system by using the 
thermodynamic curvature R of the space of thermo- 
dynamic parameters expressed through the second 
and third derivatives of the partition function or an 
appropriate thermodynamic potential. The space of 
thermodynamic states is flat if there are no interac- 
tions between particles and it curves if interactions 
are present. The thermodynamic curvature has an in- 
teresting behavior close to the critical point which 
has been shown to be singular for some well-known 
thermodynamics systems at criticality in which cor- 
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relation length is infinite (for a useful review see 
flULflEzl]). Furthermore, the thermodynamic curva- 
ture determines information about interatomic inter- 
action. Janyszek ^\ found that for an ideal classical 
gas the scalar curvature R is always zero. Janyszek 
and Mrugala [Hill showed that R is negative for Fermi 
gas and positive for Bose gas. Mirza and Moham- 
madzadeh [0, [l(J worked out the fractional statis- 
tics and also ideal q-deformed bosons and fermions 
which show the proper sign. It has also been ar- 
gued that scalar curvature is able to show that the 
fermion gas is more stable than the boson gas 111 ill . 
Other thermodynamic models, such as phase transi- 
tion properties of the van der Waals gas, have been 
investigated to find that the singular point of the ther- 
modynamic curvature is consistent with the critical 
point of the system 1112l.ll3ll. 

A standard scaling for free energy per site / in the 
critical region may be assumed as 



f(t,h) = X- 1 f(tX a ',hX 



! . ~ a h . 



(1) 



j85 



are seal- 



where, t = j8 - j8 c and a t = ^, a h - w 
ing dimensions for energy and spin operators, re- 
spectively (/3 C - designates the critical point). ^ is the 
scaling function and d is the dimension of the sys- 
tem (more details can be found in [14]. For a sec- 
ond order transition, according to the hypothesis on 



Preprint submitted to Physics Letter A 



January 15, 2013 



dimensional grounds, the curvature depends on the 
correlation volume. Since correlation length t, is the 
only physical scale in a system close to the critical 
region, one can write R ~ Combined with hyper- 
scaling, yd = 2 — a, this leads to 



R 



.a-2 



(2) 



where, a is the critical exponent related to specific 
heat. Eq. (2) is confirmed for a > [15]. It has also 
been found that this relation is not true for a < 
. [nil . In this work, we considered a case in which 
a = to find that the standard scaling would not 
work for it, either. 

The 2d Ising model in a magnetic field is a well- 
known unsolved problem in statistical mechanics. In 
1960 Fisher introduced, a remarkable solution for a 
two-dimensional super-exchange antiferromagnetic 
Ising model in a magnetic field fl 1 8M . The Fisher 
model is determined in a square lattice where there is 
an external magnetic field applied to the decorating 
spins which interact via super-exchange interactions. 
A similar super-exchange model on the kagome lat- 
tice has been considered [fl9i I20I1 . It is known that 
special cases of the kagome Ising model are soluble 
in the presence of a magnetic field 112 ll 12211 . In one 
of the special cases, this model reduces to the Fisher 
model. In this work, we calculated the scalar curva- 
ture for the Fisher model under the presence of an 
external magnetic field. 

The rest of the Letter is organized as follows. In Sec- 
tion 2, the Ising model on the kagome lattice in a 
nonzero magnetic field is reviewed. In Section 3, the 
geometrical structure and the thermodynamic curva- 
ture are introduced and the scalar curvature is calcu- 
lated for the Fisher model. The scaling behavior of 
R is also investigated. 

2. Ising model on the kagome lattice in a nonzero 
magnetic field 

Consider an Ising model shown in Fig \T\ with an 
interacting energy 

- J l Oi0 2 + J(o 2 o 3 - o 3 Oi) -iiH(o 2 + o 3 ) (3) 

around every triangle formed by the spin 0\, o 2l 03. 
Assume further that there is an external magnetic 
field, H, applied to | of the lattice sites. The reduced 




Figure 1: The kagome lattice shown in elementary cell with 
magnetic Ising spins 1,2 (solid circles) and non-magnetic Ising 
spin 3 (open circles) 



interactions and field are denoted by b = fiJ,b\ = 
/3 J \ and h = ftH, where /3 = As a result, the mag- 
netic spins 0"i and o 2 interact a super-exchange inter- 
action via the intermediate non-magnetic spin 03. 03 
is an Ising spin which is not coupled to the magnetic 
field. It has zero magnetic moment. This model does 
not admit any restrictions on any positive values of 
H, J, h, or b. The partition function of the kagome 
Ising lattice has been shown to be equivalent to that 
of a honeycomb Ising model in a zero field that may 
be solved as follows: 



Zkg 



F N Z HC 



(4) 



where, Zkg an d Zhc are partition functions of 
the kagome Ising lattice and the honeycomb Ising 
lattice, respectively, F N is an overall factor, and N is 
the number of magnetic spins of the kagome lattice 
or the number of honeycomb lattice sites. 

The per magnetic spin free energy, f KG , and the 
per site honeycomb lattice free energy, fyc, are de- 
fined as 



f KG = lim N 1 In Z KG 
f HC = lim AT 1 In Z#c 



(5) 
(6) 



Taking N — > °° limit of Eq. (4) and making use of the 
explicit expression of f H c given in Il23ll . /kg will be 



Ikg 



lnF + -ln2 
4 

2k 



+ 



\6n 2 



2n 



d6d(pZ(6,(p) (7) 



where, 

s(e,<p) 



cosh 2ri cosh 2r+\ 

sinh 2 2rcos(0 + <p) 

sinh2ri sinh2r(cos 6 +cos (p) (8) 



In the above equation, r\ and r are functions of b 
and b\. For b\ = 0, this model reduces to the Fisher 
model for the square lattice. Here, we take the Fisher 
model and consider the critical behavior of its scalar 
curvature. With b\ = 0, the partition function and the 
free energy per spin of the Fisher model in the square 
lattice are obtained to be [19] 



Ikg = 
where 

p{b,h) 



q{r) 



p(b,h)+q(r) 



-ln2 



(9) 



2 

+ ^-ln{cosh 2 /i (sinh 2 /i + cosh 2 2£)} 

2k r2n 



1 



, , , J0Jfflln{cosh 2r 
16n 2 Jo Jo 

- sinh2r (cos0 +cos<p)} (10) 



and r satisfies the equation 

_ 4r sinh 2 h + cosh 2 2b 



cosh h 



q(r) can be written as 

I r 2n i-2n 



(11) 
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d(pd 9 ln[cosh 2r 
o 

sinh2r , 

(COS0 +COS(p)}\ 



cosh 2 2r 

2n 1-2% 

/ d<pd01n{cosh 2 2r} 
Jo 

sinh2r 
cosh 2r 



(cos + cos <p) } 



This term can be further simplified to cos 6 + cos <p = 
2 cos cos If we change the variables 9\ = 



2 



,02 



6-q> 
2 



and K 
1 



2sinh2r 
cosh 2 2r ' 



we will then have 



2n 



\6it 2 jo 
+ ln[l + fccos0icos0 2 ]} 



d9i / J0 2 {ln[cosh 2 2r] 

J—n 



1 l f 2n r K 

-ln[cosh 2 2r] + -^ / dQ x \ d0 2 

2 87T Z Jo J-k 



ln[l + /ccos0i cos 62] 



1 T 1 

-ln[cOSlr2r] + - 



7t Jo 



de 



1 



ln[^(l + Vl-K 2 sin 2 )]. 



3. Thermodynamic curvature 

The geometrical structure of the statistical thermo- 
dynamic phase space was investigated in detail by 
Gibbs and later developed by Weinhold using energy 
representation and also by Ruppeiner based on en- 
tropy representation [Ql 0|. Alternatively, any ther- 
modynamic potential representation that is the Leg- 
endre transform of entropy or internal energy may 
also be used. The metric of the representation can be 
the second derivative of the thermodynamic potential 
with respect to the intensive variables. For example, 
in the Ruppeiner geometry, the extended set of r + 1 
extensive variables of the system can be denoted by 
X = (5, /V 1 , . . . , V, . . . , N r ) . The thermodynamic poten- 
tial is defined as 



4> = 4>({F''}), 



(12) 



where, $ is a Legendre transform of entropy with re- 
spect to the extensive parameter, X 1 , "r" represents 
the number of various kinds of particles and F is a 
standard intensive quantity in the entropy represen- 
tation, 

F' = (1/7/, V/7\ ...,P/T, ...n'/T) (13) 
dS (H) 



dX l 



Here, T is temperature; jl, the chemical potential; 
and P is pressure. The F l 's are a complete set of 



coordinates on the thermodynamic phase space. The 
metric in this representation is given by 



d 2 4> 
8ij ~ dF'dFJ ' 



(15) 



Another type of the Riemannian metric in the space 
of probability distribution is the Fisher-Rao metric. 
If the distribution has an exponential form such as 
Boltzmann-Gibbs distribution, the Fisher-Rao metric 
can be defined as yO, 



Sij 



d 2 \nZ 

WW 1 ' 



(16) 



where, f5' = F l /k and Z is the partition function. 
We consider a system with two thermodynamic de- 
grees of freedom. The thermodynamic surface, or the 
parameter space, is two-dimensional and the metric 
components are defined by derivative of free energy 



fij — , 



<'j 



didjf 



(17) 



where, / = InZ is the reduced free energy per site; 
and i and j represent b and h which are related to 
the inverse temperature and the external field, respec- 
tively. 

For such a metric in the two-dimensional thermody- 
namic space, the Ricci scalar may be defined by 



R 



fbb fbh fhh 
fbbb fbbh fbhh 
fbbh fbhh fhhh 



fbb fbh 
fhb fhh 



(18) 



where, fjk = didjd^f. To determine R, we need to 
obtain the second and third derivatives of free energy. 
According to the Eq. (9), the first derivatives of free 
energy are given by 



db db dr db 
dh dh dr dh 



(19) 
(20) 



Other parts of the calculations can be found in the 
Appendix. Once the components of R have been 
determined for the Fisher model, they are substituted 
in Eq. (18). The Ricci scalar is depicted as a function 
of b and h in Fig [2] From this figure it is found that at 



R(L-) o- 




b(T-') 



Figure 2: A 3d plot of R for h = 1.2 to 1.3and£ = 1 to 1.1 . 

high temperature where system is in disordered state 
(paramagnetic), the sign of R is positive (FigHJ). As 
temperature decreases R diverges to plus infinity. At 
low temperature where system is in the ordered state 
(antiferromagnetic) R is negative. When temperature 
increases R diverges to minus infinity. Ordered state 
is more stable with a lowest ground state energy. 
Curvature is singular along the critical line in the 
h — b plane which exactly corresponds to the Fisher 
expression IU9I1 as shown in Fig. |3] 



h — arccosh 



\/2-l 



arcsinh 2b] 



(21) 



Focusing on the zero field case, we have calculated R 
and its components for this special case. To find the 
critical point, the denominator of R should be zero 
and the term that makes it zero is 



cosh 2 2b + sech 2 2£ — 6 







arccosh [ 



'2+ y/2 



(22) 



where, b c is the critical point that can also be ob- 
tained from Eq. (22). The dominant components of 
R near the critical point can be written as follows 

f hb ~K(K) 
f bin fbbh i fhhh = 

f hh ~ const 

f 1 

fbbb ~ 



b-b c 
fbhh ~K(k). 



(23) 



Figure 3: A plot of h with respect to b for R 

and a plot of h = arccosh [ 
line). 



200 (gray line) 
arcsinh 2b] (black dashed 



where K(k) denotes the elliptic integral of first kind. 
In the above equations, the terms with odd numbers 
of h derivatives are equal to zero since there is ±h 
symmetry in the free energy and the other terms have 
expanded about b c to the second order. Then, R is 
given by 



R 



K(K) 

l 

b-b c 







K(K) 



const 

K(K) 





K(k) 
const 



Finally, we obtained 
1 

R 



(24) 



(25) 



K(K)(b-b c ) 

It should be noted that K( k) has a singularity at K = 1 
and can be written as 

4 



K(k) « In 



Hb-bc) 



(26) 



So, near the critical point, the scalar curvature be- 
haves as in 



R 



1 



(27) 




Figure 4: A plot of R for h = and b = 0.5 to 1 by steps 0.002. 



5=0 D24D. We conclude that ~ (b-b c )~ l = t~\ 
For h = 0, the scalar curvature R is plotted in Fig |4l 
It shows that R diverges at certain point. We care- 
fully examined this point and found that it coincides 
with the critical point. In this work, we investigated 
a soluble kagome Ising model under the presence of 
an external field in which critical exponent a is zero. 
There are standard scaling forms for a > and a < 
as R ~ t a ~ 2 and R ~ t a ~ l , respectively rfl5l ITvL [Till . 
We found that, for this model with a = in 2 di- 
mensions, the scaling behavior of R ~ 



1 is simi- 



(b-b c )ln(b-b c ) 
We can write Inx = limg^ ^ (-^ 8 ~ !)■ So, \n(b 

b c ) behaves as a power term like (b — b c ) s in which with a = (R ~ t a ~ l ). 



lar to a < of the spherical model studied in 111 711 . It 
should be noted that a = for the spherical model in 
d = 4 and van der Waals model in d = 3, the scaling 
behavior of R is ~ t a ~ 2 which is the same as a > 0. 
However in 2d-Ising model on planar random graph 
with a = — 1 and kagome lattice with a = the scal- 
ing behavior is similar as R ~ t a ~ 1 . We think that the 
scaling behavior of R is also related to dimension. 
Further studies are needed to explore this behavior. 

4. Conclusion 



The scalar curvature R, was derived for the Ising 
model on the kagome lattice in the nonzero exter- 
nal field. It was found that R has a singularity at 
critical points (Fig. 2) and that its singularity cor- 
responds to the Fisher expression depicted in Fig. 3. 
For the zero field, we derived R and its components 
close to the criticality. An expression for R was ob- 
tained for this case with a = 0. This was not what we 
expected from standard scaling of R. It has already 
been shown that standard scaling does not hold true 
for negative values of a lfl7ll . Here, we found a simi- 
lar scaling of the scalar curvature for 2-d Ising model 



(A.1) 



5. Appendix A 

Below are a few useful relations. 

dq_ 1 

db 167r(cosh4& + cosh2/z) 2 

csch[^ log(cosh 2 2bsech 2 h + tanh 2 h)] 

sech[^ log(cosh 2 2bsech 2 h + tanh 2 h)] 

{k( 1 + cosh46 + 2 cosh2/z) 2 sech 2 /z - sech 2 /z - 

16#(4sech 2 log(cosh 2 26sech 2 /z + 

tanh 2 /z)] tanh 2 log(cosh 2 2bsech 2 h + tanh 2 ft)] ) 

(cosh4Z> + cosh 2/i - 8 cosh 4 bsech 2 h smh 4h ) } 
sinh4&, 



(A.2) 



^ _ i 

dh 167r(cosh4& + cosh2/z) 2 
csch[^ log(cosh 2 2bsech 2 h + tanh 2 h)] 

sech[- log(cosh 2 2Z>sech 2 /j + tanh 2 h)\ 

{k{\ +cosh4£ + 2cosh2/i) 2 sech 2 /z- 

16i^(4sech 2 [^log(cosh 2 2bsech 2 h 

+ tanh 2 /z)] ) (cosh4l? + cosh 2h - 

8 cosh 4 &sech 2 /z sinh 4 *) } sinh 2 2Z?tanhr(A.3) 
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dp 
db 



cosh 2b sinh 2b 
cosh 2 2b + sinh 2 h 
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